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ABSTRACT 



We report observations from the Spitzer Space Telescope (SST) regarding the 
frequency of 24 /im excess emission toward sun-like stars. Our unbiased sample 
is comprised of 309 stars with masses 0.7-2.2 M and ages from <3 Myr to >3 
Gyr that lack excess emission at wavelengths <8 /xm. We identify 30 stars that 
exhibit clear evidence of excess emission from the observed 24/8/im flux ratio. 
The implied 24 /xm excesses of these candidate debris disk systems range from 
13 % (the minimum detectable) to more than 100% compared to the expected 
photospheric emission. The frequency of systems with evidence for dust debris 
emitting at 24 /j,m ranges from 8.5-19 % at ages <300 Myr to < 4 % for older 
stars. The results suggest that many, perhaps most, sun-like stars might form 
terrestrial planets. 

Subject headings: planetary systems: formation; infrared: stars; stars: circum- 
stellar matter 
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Introduction 



Are planetary systems like our own common or rare in the Milky Way galaxy? The 
answer depends on what aspect of our planetary system one is investigating. Gas and 
dust rich circumste llar disks appear to be a common outcome of the star formation process 
( jStrom et al.lll993l ). Gas giant planets within 5 A U (presumably formed from these disks) 
surround >6% of sun- like stars (jMarcv et al.l 120051 ) while the detection of terrestrial planets 
is still in its infancy (IBeaulieu et al.l 120061). Al t hough debates concerning t heories of giant 



planet formation continue (IDurisen et al 



2007: Lissauer fc Stevenson 20071), there is some 



consensus regarding the formation of terrestrial planets (e.g. iNagasawa et al.ll2007l ). Starting 
with a swarm of 1 km-sized planetesimals, orderly growth of larger bodies proceeds rapidly 
(< 1 Myr) out to at least 2 AU. When the gravitational cross section greatly exceeds the 
geometrical cross section of the largest objects, growth transitions from orderly to oligarchic 
with the biggest bodies growing fastest in a runaway process. The final stage, chaotic 
growth, is characterized by high velocity collisions between the few remaining large bodies 
in the system. Remaining challenges incl ude the formation of km-sized planetesimals in 
the face of gas drag on meter-sized bodie s (lWeidenschillindll977l ) and Type I migration of 
lunar-mass objects in a remnant gas disk (INelsonl 120051 ). 



Yet there are few observational tests of this developed theory. The physical character- 
istics of the terrestrial planets, their satellites, and the asteroid belt provid e constraints on 
the formation of our solar system (IBottke et al.l 120051 ; lO'Brien et al.l 120061 ) . Observations 
of circumstellar dust debris surrounding su n-like stars can be used to trace the presence of 
planetesimal belts of larger parent bodies (IMeyer et al.l 120071 ) and thus constrain theories 
of planet formation. Far-IR observations at 70 /zm with the Spitzer Space Telescope sug- 
gest that 10-15% o f sun-like stars poss ess cool outer dust disks that are massive analogs 
of the Kuiper Belt ( iBryden et al.l 120061 ) . Yet few mature stars exhi bit mid-infrared excess 
indicative of terrestrial temperature material (IBeichman et al.l 120051 ) . Recent work with the 
Spitze r Space Telescope has begun to assess th e frequency of this emission around younger 
stars (jChen et al.l 120051 ; iHernandez et al.ll2006l ). In this contribution we use Spitzer data 
to investigate the frequency of mid-IR excess emission, which may originate from 1-10 AU, 
observed toward sun-like stars over a wide range of ages spanning the epoch of terrestrial 
planet formation in our Solar System. 



2. Observations 



Observations were obtained a s part of the Forma tion and Evolution of Planetary Systems 
(FEPS) Legacy Science Program (IMeyer et al.l 120061 ). Our sample consists of 328 "sun-like" 
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stars with spectral types F5-K3 and masses ranging from 0.7-2.2 M Q (though strongly peaked 
at 1.0 M ). The sample was constructed so that roughly equal numbers of stars were selected 
in logarithmically spaced age-bins from 3 Myr to 3 Gyr (each bin spanning a factor of x3 in 
age). Stars < 100 Myr were largely drawn from young stellar populations within the Local 
Association, often members of OB and T Associations. Ages for these stars were estimated 
from pre-main seq uence evolutionary tr acks, as well as kinematic association with groups 
of known age (e.g. iMamajek et al.ll2002l ). Older main sequence stars were selected from a 
volume-limited sample of stars taken from the HIPPARCOS catalog. Ages for these stars 
were estimated from c alcium H &: K e mission-line indices which trace stellar activity levels 
using the calibration of lDonahud (119981 ) . Errors in age for both the young stellar populations 
as well as the main sequence stars are estimated to be < 50 %, though uncertainties in the 
absolute calibration of these ages are not well understood. A fraction of our sample in the 
age range from 30 Myr to 1 Gyr were s elected to be members of open clust ers with well 



determ ined ages (e.g. IStauffer et al.ll2005l ). Our sample selection is described in lMeyer et al. 



(120061 ) and details concerning the age estimates for each star are given in Hillenbrand et al. 
(in prep). We fit Kurucz model atmospheres to B/V (Tycho) photometry from HIPPARCOS 
and JHK5 photometry from 2MASS. We assumed solar metalicity and surface gravities 
estimated from the position of each star in the H R diagram, performing a non-linear least 
squares fit for T e ff and solid angle. For stars within 75 pc, we assumed A v =0 while for more 
distant targets A v was a free parameter. 

All sources were detected at 8 fim and 24 /zm with SNR > 30 using the SST. Photometry 



at 8 /im was derived from sub-array observations with the IRAC instrument (IFazio et al. 



20041 ) . We began our analysis with data processed through the S13 pipeline. Cosmic ray 
rejection was implemented and corrections for spatially-dependent pixel area and filter re- 
sponse variations were applied. Aperture photometry was derived from from a 3.7" radius 
aperture (with sky anuli ranging from 12.2-24.4") usi ng a modified vers ion of IDLPHOT 
and placed on the standard flux scale recommended by lReach et al.l ( 120051 ). Systematic cal- 
ibration errors are estimated to be < 2 % in each band. Random photometric uncertainties 
were estimated from the repeatability of 64 observations obtained at each dither position for 
each source resulting in minimum uncertainties of 1% at 8/im. Photometry at 24 /mi was 
derived from eithe r 28 or 56 exposur es with integration times of 3 or 10 seconds each with the 
MIPS instrument (IRieke et al.ll2004l ). We began with S13 pipeline data and photometry was 
derived using the MOPEX software. Fluxes were estimated from a PSF-fitting algorithm 
and placed o n the flux scale recomm ended by the SSC. Calibration errors are thought to 
be < 4% (cf. lEngelbracht et al.l 120071 ). The minimum random uncertainties in the 24 /im 
photometry are 1 



70. 
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3. Analysis 

Of the parent sample of 328, 14 sta rs were selected for our initial IRS search for remnant 



circumstellar gas (IPascucci et al.ll2006l ) on the basis of previously detected dust disk signa- 
tures. Our analysis starts with the unbiased FEPS sample of 314 (= 328 - 14) stars spanning 
a range of ages from <3 Myr to >3 Gyr. We use the 24/8/im flux ratio to search for stars 
that exhibit excess emission. These data are plotted in Figure 1 as a function of 8/im flux. 
Note that the brighter sources tend to be nearby (older) field stars, while fainter targets tend 
to be more distant (and younger) sources. Several systems exhibit flux ratios indicative of a 
24 /im excess. The expected photospheric ratio (24/8/im) in this diagram is approximately 
0.116. In order to define an empirical "blue envelope" of stars without excess, we employed 
a sigma-clipping algorithm to the distribution of 24/8/im flux ratios. Initially, the mean 
flux ratio was computed to be 0.190 with a = 0.82, and two outliers with ratios beyond 
3a (indicating obvious excess). Removing these two outliers, we recompute the mean and 
sigma, resulting in identification of three additional sources with smaller excesses. Repeating 
this process a total of seven times, the values converge with a mean of 0.117 and a = 0.005 
as shown in Figure 1, consistent with model predictions. We identify 35 (positive) outliers 
which we attribute to excess emission in the 24/im b and. Of these, only five exhibit excess 



emission at < 8/im and were previously identified in ISilverstone et al.l (120061 ) as optically- 
thick primordial gas rich disks. Because we are interested in understanding the transition to 
debris disks at 24 /im we remove these five from the sample of excess stars under considera- 
tion (314 — 5 = 309). Assuming that the color excesses observed are due to excess emission 
at 24 /im, our 3a detection limit of (0.117 + 3 x 0.005 = 0.132) corresponds to a 13 % excess 
at 24 /im compared to the expected photospheric emission (0.132/0.117 — 1). The largest 
inferred 24 /im excess is just over 100 % compared to the photosphere (0.245/0.117 — 1). The 
30 sources with detectable excess from our sample of 309 are listed in Table 1 as a function 
of age. According to the Shapiro-Wilk test, the distribution of 24/8/im flux ratios for the 
279 sources without excess is not gaussian (P < 0.01 %). We tested to see whether the mean 
ratio of 24/8 /im emission was a function of source brightness. For the 140 targets with 8/im 
> 128 mJy the mean was 0.1137 with a = 0.0036, while for the 139 fainter targets the mean 
was 0.1199 with a = 0.0047. This small offset could be due to uncertainties in flux calibra- 
tion as a function of integration time (Carpenter et al. 2007; Engelbracht et al. 2007)0. As a 
result, the errors quoted on the reported excesses include the random errors in the 24/8/im 
ratio (typically 1-2 %), as well as the dispersion in our estimate of the photospheric color 
(4.3 %) rather than the error in the mean, added in quadrature. 



1 Adopting these offsets in the mean flux ratio (and associated 3 a limits) would result in identification 
of one new excess object (HD 43989, 30-100 Myr old) and removal of [PZ99] 161618.0-233947 from Table 1. 
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In Figure 2, we present the fraction of stars exhibiting 24 /zm excess emission in our 
sample as a function of age. Eac h bin sp a ns a f actor of 3 in age. The errors in the ordinate 
are Poisson, computed following Gehrelsl ( 19861 ) with excess fractions as follows: (5/30) for 
stars 3-10 Myr, (9/48) for 10-30 Myr, (5/59) for 30-100 Myr, (9/62) for 100-300 Myr, (2/53) 
300-1000 Myr, and (0/57) for stars 1-3 Gyr old. The KS test suggests that the distributions 
of 24/8/im flux ratios (Fig. 1) for the sample < 300 Myr (N = 199) and those > 300 Myr 
(N = 110) are inconsistent with having been drawn from the same parent population (P < 
10~ 10 ). We note that the errors in age quoted above act to diffuse sources to younger as well 
as older ages. Because there are more excess objects in younger bins, errors in age tend to 
increase the excess fractions at older ages. As a result, the abrupt drop in the excess fraction 
at 300 Myr may be even more dramatic than detected here. 



Discussion 



We associate the observed 24 /im excess with dust debris generated through collisions 
of planetesi mals. One of our excess stars identified in Table I (HD 12039) was studied 
in detail by iHines et al.l (120061 ). Models of this debris disk (with fractional 24 /im excess 
0.151/0.117 - 1 = 0.29) suggested a dust mass of ~ 10~ 5 M e located between 4-6 AU. 
The magnitude of all our detected 24 /im excesses are within a a factor of x 3 (relative to 
the photosphere) compared to HD 12 039. Results similar to ours, have been reported for 
samples of FGK stars in open clusters (IGorlova et al.ll2006l ; ISiegler et al.l 120071 ). Our sample 
is comprised of 60 open cluster stars with discrete ages of 55 (5 members of IC 2602), 90 (13 
members of a Per), 110 (20 Pleiades), and 600 Myr (22 Hyades), as well 249 field stars. We 
have analyzed the statistics for sub-samples where they overlap. While the excess fractions 
for open clusters with ages 30-100 Myr and 100-300 Myr are greater than the comparable 
field star samples(3/18 vs. 2/41 and 5/20 vs. 4/42 respectively), the results are formally 
consistent with each other. This suggests that there is no strong dependence of debris disk 
evolution on star-forming environment, though larger samples could reveal a difference. 

A key question is whether stars observed to have excess at 24 /im in one age bin are the 
same cohort of stars with 24 /im excess in another. In other words, do the same 10-20 % 
of sun-like stars with excess evolve from one age bin to the next with a constant fraction; 
or are they distinct groups of stars, that persist in the observed state for a short time? Our 
observations trace excess emission from 21.7-26.4 /im toward stars lacking excess emission < 



10/qr i (corresponding to a lack of dust generating planetesimals inside 1 AU (jSilverstone et al. 



20061 )). Assuming blackbody emission from large grains implies dust at radii from ~ 4-7 AU. 
Maximum dust production during the evolution of a planetesimal swarm is thought to occur 
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between runa way and chaotic growth when t he largest planetesimals reach ~ 2000 km at a 



given radius (IKenyon fc Bromleyl I2004J . 120061 ) . The timescale for this goes as r 



disk 



where a is the orbital radius and Udisk is the mass surface density of so lids in the disk 
(IGoldreich et al.l (120041 )). Assuming o^isfc ~ c a~ p , and that < p < 1 ( IKitamura et al. 



( 20021 )). a range of x2 in radius corresponds to x3 — 6 in time. Thus the 



emission we 



observe might not persist over timescales much larger than our age bins q Perhaps many 
stars go through this phase of 24 /xm excess, but a t diff erent times. A range of x 100 in 
initial mass surface density (jAndrews Williams! (120051 ) ) could translate into a range of 
x 100 in evolutionary timescales. If so, one might expect smaller excesses at later times 
(produced by lower mass disks). In comparing the mean detected excess for stars 3-30 Myr 
old (0.359 with a = 0.199) with that for stars 30-300 Myr old (0.345 with a = 0.116), we 
find no evidence that this is the case (though the samples are dominated by stars lacking 
detectable excess). Nevertheless, one might consider summing the fractions of stars with 
24 (Ltm excess between 3-300 Myr, resulting in an overall fraction of stars with evidence for 
terrestrial planet formation greater than 60 %! Averaging the results over factors of ten in 
age results in excess fractions of 18, 12, and 2 % at ages 3-30 Myr, 30-300 Myr, and 0.3-3 
Gyr, implying at least 32 % of sun-like stars exhibit evidence for terrestrial planet formation 
(provided that the epoch of 24 /im excess emission lasts < xlO in age). We note that in this 



scenar io, the planets formed later from lower mass disks will be smaller (IKenyon &: Bromley 



Results to date suggest that: a) primordial disks between 0.3-3 AU dissipa te or agglom- 
erate into larger bodies on timescales comparable to the cessation of accre tion Jrlaisch et al. 



2001); and b) few stars harbor optically-thin inner disks between 3-30 Myr (ISilverstone et al. 
20061 ). Based on the oretical con si derat ions, we expect that planetesimals belts evolved 



rapidly within 3 AU. iRieke et al.l (120051 ) explore the evolution of 24 /^m excess emission 
around a sample of A stars observed with SST and IRAS. They deduce a characteristic 
timescale of 150 Myr for strong excesses to decay. However, care must be taken in compar- 
ing these results to ours as: 1) the dust masses detected by Rieke et al. are likely larger 
than the dust masses detected here; and 2) similar temperature dust traces distinct radii for 
stars of different luminosity. In general, the fractional 24 fiia excesses around A stars are 
larger than around G stars. The observed duration of the excess phase for both samples is 
longer than expected if the emission results solely from dust production well inside 10 AU. 
An important caveat to our results is that we have not assessed whether the 24 /im excesses 



2 While the published Kenyon and Bromley models predict the duration of 24^m excess emission > x 3 — 6 
in time, they also predict hot dust at smaller radii covering a wider range of radii than our observations 
imply. 
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we have detected are tracing warm dust in the terrestrial planet zone, or the Wien-side of 
the Planck function from cooler dust. If, at 24 /im, we are seeing cooler dust generated 
at radii beyond 10 AU, we would expect to observe it at later times. We also note that 
the maximum excess ratios predicted by Kenyon and Bromley are larger than the excesses 
detected here. 



Wyatt et al.l (120071 ) provide steady-state models of warm dust production around sun- 
like stars. On the basis of comparing the observed IR luminosity of several systems to these 
models as a function of age, they conclude that most (5/8) of the stars exhibiting evidence 
for warm dust must be in a transient state of evolution and not participating in a steady- 
state collisional cascade. However, the three systems with ages < 300 Myr (one of which is 
HD 12039) can be explained with equilibrium models. Global analysis of the SED for all 
sources identified here will be required to directly compare the evolutionary state of these 
systems with the models (Carpenter et al., in preparation). 



More work is needed to define the transition from primordial to debris disk (e.g. lPadgett et al 



20061 ). Given the short time expected for collisional evolution of inner planetesimals belts 
(< 3 Myr), and the 1-10 Myr lifetime of primordial disks, it may be difficult to detect the 
onset of collisional evolution. We suggest that SST observations at 24 /im can be interpreted 
as evidence for terrestrial planet formation occuring around many (19-32 %), if not most 
(62 %), sun-like stars. This range is higher tha n the observed frequency of gas giant planets 
(6.6-12 % within 5-20 AU iMarcy et al.l (120051 )) but comparable to the inference that cool 
dust debris beyond 10 AU might be very common iBryden et al.l (120061 ). Radial velocity mon- 
itoring of low mass stars, micro-lensing surveys, as well as transit surveys such as COROT 
and Kepler, will provide critical tests of our interpretation. 

We thank all members of the FEPS, IRAC, MIPS, and SSC teams for their efforts, 
as well as Scott Kenyon, Nick Siegler, and George Rieke for valuable discussions, and an 
anonymous referee for helpful suggestions. This work was based on observations made with 
the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA. Support for this work was provided 
by NASA through an award issued by JPL/Caltech. 
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Table I. Systems with MIPS-24 /im Excess 



Source 


Dist 




T 


logf'L^ /Lin ) 


fodurr, (cXCCSs) /fodiirn I'DllOt) 




(J 




(pc) 


(yr) 


(K) 


(dcx) 








1RXS J051111. 1+281353 


199 


6.5-7 


5270 


0.71 


0.239 





.045 


RX J1600.6-2159 


161 


6.5-7 


5330 


0.27 


0.190 





.045 


[PZ99] J161459.2-275023 


114 


6.5-7 


5500 


-0.14 


0.628 





.047 


[PZ99] J155847.8-175800 


161 


6.5-7 


4660 


0.20 


0.479 





.046 


[PZ99] J161618.0-233947 


161 


6.5-7 


5250 


0.22 


0.141 





.045 


HD 22179 


68 


7-7.5 


5990 


0.02 


0.336 





.045 


HD 116099 


140 


7-7.5 


6010 


0.19 


0.188 





.045 


HD 141943 


67 


7-7.5 


5810 


0.43 


0.210 





.045 


HD 281691 


73 


7-7.5 


5140 


-0.42 


0.156 





.045 


MML 8 


108 


7-7.5 


5810 


0.15 


0.737 





.047 


MML 17 


124 


7-7.5 


6000 


0.43 


0.612 





.046 


MML 28 


108 


7-7.5 


5000 


-0.35 


0.413 





.046 


MML 36 


98 


7-7.5 


5270 


0.03 


0.541 





.046 


MML 43 


132 


7-7.5 


5410 


0.06 


0.154 





.045 


HD 377 


39 


7.5-8 


5850 


0.09 


0.332 





.045 


HD 12039 


42 


7.5-8 


5690 


-0.05 


0.287 





.045 


HE 750 


176 


7.5-8 


6360 


0.28 


0.197 





.046 


HE 848 


176 


7.5-8 


6310 


0.47 


0.537 





.047 


W79 


152 


7.5-8 


5380 


-0.29 


0.242 





.046 


HD 19668 


40 


8-8.5 


5420 


-0.23 


0.192 





.045 


HD 61005 


35 


8-8.5 


5460 


-0.26 


1.096 





.049 


HD 72687 


46 


8-8.5 




-0.05 


0.208 





.046 


HD 107146 


28 


8-8.5 


5860 


0.02 


0.329 





.045 


HII 152 


133 


8-8.5 


5700 


-0.10 


0.366 





.046 


HII 250 


133 


8-8.5 


5770 


-0.04 


0.146 





.046 


HII 514 


133 


8-8.5 


5720 


0.04 


0.250 





.046 


HII 1101 


133 


8-8.5 


6070 


0.08 


0.545 





.046 


HII 1200 


133 


8-8.5 


6210 


0.35 


0.170 





.045 


HD 85301 


32 


8.5-9 


5600 


-0.14 


0.343 





.046 


HD 219498 


62 


8.5-9 


5670 


-0.07 


0.277 





.045 
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Fig. 1. — 24/8/im flux ratio plotted as a function of 8/im flux for 314 stars drawn from 
the unbiased FEPS sample. The sample mean (solid line) and 3a limits (dashed lines) as 
described in the text are shown. 
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log(Age yrs) 

Fig. 2. — The fraction of stars in the sample with detectable 24/zm excess plotted as a 
function of age. 



